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The  Accuracy  of  Remapping  Irregularly  Spaced  Velodty  Data 
onto  a  Regular  Grid  and  the  Computation  of  Vorhcity 

R  K.  Cohn  and  M.M.  Koochesfabani 

Department  of  Mechanical  Engineering,  Michigan  State  University,  East  Lansing,  MI 
The  velocity  data  obtained  from  Molecular  Tagging  Velocimetry  (MTV)  are  typically  located 

by  photons.  The  velocity  vector  is  determined  from  the  displacement  0  *  J  This 

a  nulsed  laser  which  are  imaged  at  two  successive  tunes  within  the  lifetime  of  A*  tracer.  this 
technique  may  be  viewed  as  the  molecular  counteipart  of  PIV.  To  take  advantage  ofsta^d^  ^ 
Drocessing'teclmiques,  the  MTV  data  need  to  be  remapped  onto  a  regular  gnd  with  a  unrform 
spacing  li  this  work  we  examine  the  accuracy  and  noise  issues  related  to  the  use  of  various i  low 
order  polynomial  least-square  fits  for  remapping  and  the  subsequent  computation  of  vorticityfrom 
these  data.  The  information  obtained  has  relevance  to  PIV  data  processing  as  we  .  y 

Spedding  and  Rignot  (1993),  the  best  estimate  of  the  location ^of  the velocity  vector ^^uired  tluou^ 
the  use  of  tracer  techniques,  such  as  PIV,  is  at  the  midpoint  of  the  displacement  vector.  Thus,  nl 

special  care  is  taken,  PIV  data  are  also  initially  obtained  on  an  irregular  gnd. 

In  the  past,  various  methods  have  been  used  for  remapping  randomly  spaced  velocity  data  onto 
a  regular  grid.  Among  them  are  an  inverse  distance  approach  and  a  “global  basis  funcbon 
by  Spedding  and  Rignot  (1993).  In  this  study,  we  consider  the  use  of2  ,3  ,or4  order  pol>Tio 
to  remap  the  velocity  field  by  performing  a  local  least-squares  fit  of  die  nreguhtfly  spaced^ 
within  region  of  radius  R.  This  choice  was  selected  based  on  the  work  of  Agui  and  Jimenez  (198  /), 
who  report  that  low  order  polynomial  fits  and  “kriging"  techniques  produce  the ^most  aeon-are 
representation  of  the  actual  velocity  field.  However,  no  quantitative  information  on  the  performance 
of  these  methods  is  presented*  Four  approaches  are  assessed  here  for  computing  the  ou  -o  p  ane 
vorticity  field  from  the  in-plane  velocity  measurements:  direct  differentiation  of  the  polynomial  fits 
used  in  the  remapping  process,  1st  and  2nd  order  finite  difference  techniques  (2  mid  4  order 
accurate,  respectively)  and  an  8-point  circulation  method  on  the  regular  data.  It  should  be  noted  that 
the  8-point  circulation  method  is  identical  to  the  1“  order  finite  difference  method  after  an 

appropriate  choice  of  smoothing.  _ 

Several  authors  have  examined  the  accuracy  of  the  various  means  to  compute  vorticity  from 
regularly  spaced  data.  Spedding  and  Rignot  (1993)  used  a  1-  order  finite  difference  technique  for 
the  inverse  distance  method  or  directly  differentiating  the  global  basis  function.  Results  indicated 
that  the  global  basis  function  produced  generally  superior  results,  however,  results  were  found  to  be 
highly  dependent  upon  the  ratio,  of  a  characteristic  length  scale,  of  the  flow,  L,  to  the  mean  sp^ng 
between  measurements,  6.  Abrahamson  and  Lonnes  (1995)  found  that  the  circu  ation  me  o 
resulted  in  slightly  more  accurate  vorticity  results  than  using  a  least-squares  fit  to  a  model  velocity 
field.  Luff  et  al (1999)  compared  the  1st  and  2Dd  order  finite  difference  methods  and  an  8-pomt 
circulation  method  in  the  calculation  of  vorticity  in  the  presence  of  both  noise  and  missing  data 
points.  In  terms  of  only  the  computed  vorticity  RMS,  the  1st  order  finite  difference  technique 

produced  the  best  results.  .  , 

One  short-coming  of  the  above  mentioned  studies  is  that  only  the  random  component  ot  the 

error  field  is  examined.  Fouras  and  Soria  (1998)  found  that  the  error  in  the  vorticity  field  could  be 
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better  represented  if  it  is  divided  into  two 

the  fit  This  wotk  was  based  entirely  on  regularly  sampled  velocity  data,  issues  conn 

™3o 

denendtag  on  the  criterion  used  to  assess  the  error.  In  our  work  we  directly  compare  several  of  d» 
different  vorticity  calculation  methods  which  were  determined  in  the  previous  studies  to  produce 
SSStaSSKn.  we  include  the  effect  of  remapping  of  the  velocity  field  on  the  accuracy  o 
vorticitv  estimation.  The  differentiation  of  the  polynomial  fit,  used  in  remapping,  as  a 
Ista  Srvortci.y  is  also  considered.  A  simulation  of  an  Oseen  vortex,  utilized  by  the  previous 
“  ae  ‘  if  also  employed  here  as  the  basis  for  comparison.  The  effect  of  uncertainty  m  the  ve  oclty 
measurements  in'sirnulated  by  adding  a  mndorn  amount  of  noise  to  each  velocity  componem. 

Our  simulations  show  that  the  accuracy  of  the  remapping  process  depends  on  both  I/d ' “d 
„onn“siae  of  the  region  used  in  the  least-squares  fit,  M.  The  effect  of  am-V  the  value 
of  S/8  on  the  velocity  and  vorticity  fields  is  to  decrease  the  random  error,  but  increase  the  mean  bias 
error.  Generally,  the  decrease  in  random  error  is  smaller  thjtn  the .increase  in the bras 
the  different  polynomial  orders  tested  produced  accurate  results  in  the  remappi  g  ty 

for  SSTSS  of  Ui  and  J»  .  For  example,  with  1/5  >  4.5  both  the  velocity  mean  bias  etror 

and  random  error  are  less  than  1%  of  the  peak  velocity.  directly 

Results  show  that  the  most  accurate  vorticity  results  are  ■**“**&  by  ^er^^Uy 

differentiating  the  3rd  order  polynomial  fit  to  the  origin  megu  ar  a  a  because  thev  either 
difference  technique.  Note  that  this  conclusion  is  different  from  previous  studies  because  they  either 

did  not  consider  the  2°“  order  finite  difference  approach,  or  the  performance  was  based  ody  on  the 
random  error.  For  i/5  >  4.5,  we  find  the  mean  bias  error  and  random  error  m  vorticity  can  be  less 
than  3%  and  2%,  respectively,  of  the  peak  vorticity.  A  decrease  in  Uh  causes  a  large grease  in  e 
vorticity  bias  error,  error  values  higher  than  18%  of  the  peak  vorticity  are  found  for  L! 6  -  2.5. 
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